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FOREWORD 


This is the final report for contract NAS8-38856 from NASA Marshall Space 
Flight Center. The study contract title was "Structural Damage Prediction and 
Analysis for Hypervelocity Impacts." The Technical Monitors were Dr. Joel 
Williamsen, Greg Olsen, and Jennifer Robinson. The work was performed 
between October, 1990 and October, 1995. The bulk of the analysis was performed 
through 1992, while continued testing and debris cloud models continued through 
the end of the contract in October , 1995 
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INTRODUCTION 


It is necessary to integrate a wide variety of technical disciplines to provide an 
analysis of structural damage to a spacecraft due to hypervelocity impact. There 
are many uncertainties, and more detailed investigation is warranted, in each 
technical discipline. However, a total picture of the debris and meteoroid hazard 
is required to support manned spaceflight in general, and the international Space 
Station in particular. In the performance of this contract, besides producing a 
handbook, research and development was conducted in several different areas. 

The contract was broken into six separate tasks. Each task objectives and 
accomplishments will be reviewed in the following sections. The Handbook and 
separate task reports are contained as attachments to the final report. The final 
section summarizes all of the recommendations coming out of this study. 

The analyses and comments are general design guidelines and not necessarily 
applicable to final Space Station designs since several configuration and detailed 
design changes were being made during the course of this contract. Rather, the 
analyses and comments may indicate either a point-in-time concept analysis, 
available test data, or desirable protection goals, not hindered by the design and 
operation constraints faced by Space Station designers. 


TASK 1 - PROBABILITY ANALYSIS 

1.1 Objective 

The objectives of this task were to develop a computer program (and associated 
user’s manual) to analyze the probability of penetration and structural damage, 
and to perform analyses for the original Space Station Freedom design. 

1.2 Results 

The analyses and comments are general design guidelines and not necessarily 
applicable to final Space Station designs since several configuration and detailed 
design changes were being made during the course of this contract. 

A series of computer programs, with the collective name SD_SURF (for Space 
Debris SURFaces) were developed under this contract. The computer programs 
and user's manuals are Attachment 1. An AIAA paper describing SD_SURF is 
included in attachment 2. The SDJ5URF codes calculate which impact velocities 
and obliquities most influence the probability of no penetration (PNP) or critical 
damage (PNCD) of a Space Station module. 

SD_SURF contains FORTRAN programs which operate with the BUMPERII 
computer code, developed by Boeing for NASA-Marshall Space Flight Center. 
SD_SURF was maintained through two major updates to BUMPERII. 

SD_SURF was compiled and run on a VAX computer and also on Apple 
Macintosh computers. Source code is supplied for each as well as compiled 
applications for the Apple Macintosh. BUMPERII was also compiled to run on 
the Macintosh computer, with minor modifications from the VAX version. 
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Several BUMPERII problems and potential improvements were identified and 
sent to the responsible engineers at NASA-MSFC. These reports are given in 
Attachment 3. Additional recommendations to improve the computer programs 
are included at the end of this report. 

SDJSURF for EXCEL contains Microsoft EXCEL spreadsheet templates and 
macros which may operate either with the output of the FORTRAN programs or 
as standalone applications. Despite its slower speed, SD_SURF for EXCEL has 
several advantages over the FORTRAN version: 

• It is easier to customize ballistic limits, or the space debris flux. 

• There is easy access to two and three dimensional graphing. 

• Simple geometries may be built in SDJSURF without the need to use 
PATRAN or IDEAS (SuperTab). 

1.3 Conclusions 

Figure 1 shows probability of penetration as a function of impact velocity and 
obliquity. The velocities greater than 9 to 14 km/s and obliquities from 40 to 65 
degrees, have the biggest influence on the PNP or PNCD of a Space Station 
module. The SD_SURF calculation was for an early Space Station Freedom (SSF) 
shield design (not the current international Space Station module design): 

• bumper: 0.127 cm 6061 aluminum 

• standoff of the bumper from the rear wall: 10.16 cm 

• Intermediate catcher: MLI with one layer of beta cloth 

• pressure vessel: 0.318 cm 2219-T87 aluminum 

Early estimates of the projectile diameter (as a function of velocity and 
obliquity) to cause critical damage to SSF, indicated that if penetration should 
occur, there was a 10 to 25 percent chance that complete rupture would occur. 

The analysis depends on final design, impact conditions, and how much 
momentum can actually be transferred to the pressure vessel wall. If penetration 
is easier than the analysis prediction at high velocity and obliquities (if less 
momentum is deposited in the rear wall), then rupture could not occur. 

The SSF wall design, with MLI half the distance to the bumper was better than 
previous designs with MLI closer to the rear wall, since this will minimize 
momentum absorbed by the rear wall. Deleting the beta cloth, or moving the MLI 
closer to the bumper would further reduce the risk of rupture, but with an 
associated reduction in PNP. However, to tell if there is a significant increase in 
PNCD, more information on high velocity impacts is required. Until then, it 
should be assumed (perhaps conservatively) that melting and vaporization of the 
impacting particle may lead to a burst type failure at high velocities. 

To reduce the PNP a heavier bumper and intermediate shield and rear wall 
(0.46 cm) were used in the final design of the U.S. laboratory and habitation 
modules for the international Space Station. This total shield redesign 
significantly reduced the chance of critical damage. The heavier shield 
(including the rear wall) reduced the initial crack length that can be generated by 
a penetration, and it increased the critical flaw size for rupture. The probability of 
critical damage became more dependent on pentration effects. 
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Flux (which penetrates) • Exposure Area * Time 

= !n(1/PNP) 


BUMPERII 

* Thick. Bumper [inch] 0.050 

* Thick. RW [inch] 0.1 25 

* Spacing [inch] 4 

GEOMETRY 


• MB17-AII.UNI 

• One Module 

-12 degrees pitch 

• Threats in GEOM 45 

• Material IDs 1,2 

• 0.5 km/s Intervals 

ENVIRONMENT 

• SSP 30425 

•Year 1995 

• Inclination 285" 

• Altitude [km] 388 

l • Mass growth rate 0.05 

• Flux factor 4 

• Solar Flux Fact 100 

• Time [yr] 1 
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Flux of Penetrating Debris vs. Impact Velocity and Obliquity 


Fig. 1. Probability of penetration as a function of velocity and obliquity. 

SD-SURF calculation for an early Space Station Freedom design. 
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TASK 2 - FRAGMENT PENETRATION / DEBRIS CLOUD MODELING 
2.1 Directives 

The objectives of this task were to model the debris cloud formed by a 
hypervelocity impact , and to define the ballistic limit surface for fragment 
penetration. 

2J2 Results 

Attachment 4 contains a summary debris cloud radiographic measurements 
by Mr. Andrew Piekutowski of the University of Dayton Research Institute. This 
addresses both the largest residual fragment as well as the shape and 
momentum distribution of the debris cloud. Two papers prepared by Mr. 
Piekutowski are contained in Attachments 5 and 6. The first, described the 
largest fragment produced by a hypervelocity impact, was presented at the 1992 
AIAA Space Programs and Technologies Conference in Huntsville, Alabama. 

The second has been submitted for review at the 1992 Hypervelocity Impact 
Symposium. Further review of the fragment penetration ballistic limit, and 
additional test data is contained in the Handbook. 

2.3 Conclusions 

The descriptions and analysis of post impact debris cloud radiographs provides 
a valuable benchmark for hydrocode analysis, as well as a test derived cloud 
model for momentum analysis of rear wall deformation. 

Three fragment types are of concern. For spherical projectile, either the 
largest remaining fragment or a general dispersion of uniformly sized fragments 
can produce penetrations. For cylindrical projectiles, a fragment of spall from 
the bumper can be the most lethal. This bit of spall comes off normal to the 
impacting projectile flat, and it is faster than the remainder of the debris cloud. 

The penetration resistance of shields with intermediate shields (including 
MLI with beta cloth) requires empirical determination. Projectile shape can have 
a significant influence on the ballistic limit. However, the SD_SURF computer 
code showed that the fragment penetration mechanism has very low influence on 
the overall PNP for current SSF designs. 

The mass and velocity distribution that impacts the rear wall is strongly 
influenced by the presence of intermediate shields. The effects are modeled 
analytically, but there is no experimental verification. 
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TASK 3 - STRUCTURAL ANALYSIS 

3.1 Objectives 

The objective of this task was to analyze structural deformation and fracture by 
an impulsive load generated by a hypervelocity impact. This includes 
deformation, crack growth and residual strength assessment. 

3.2 Results 

Rear wall thickness and permanent displacement was measured on several 
test articles that had sufficient Kevlar doth intermediate shields to prevent 
fragment penetration. These results are presented in die Handbook for 
Structural Damage Prediction and Analysis for Hypervelocity Impact 
(Attachment 7). 

Mr. David Grove of UDRI performed EPIC analyses of rear wall deformation 
due to an initial Gausian momentum intensity distribution (Wilkinson's 
assumption). He curve fit these and previous STEALTH analysis results to give a 
dimensionless parameter to predict failure. This work is summarized in the 
Handbook and included in detail in an Appendix to the Handboo. 

Mr. William Roberts performed HULL and NASTRAN analyses at the 
Michoud Assembly Facility for comparison to the measured deformation. The 
HULL analyses gave the correct qualitative shape, although more analyses would 
be required to completely match the correct momentum (with rebound), 
momentum distribution, and material properties. The NASTRAN 3D sheet 
elements appeared to accentuate the thinning and deformation at the center of the 
specimen. This was true for 10.16 and 30.48 cm radius pie slices. 

Finite element models of a SSF size cylinder with a longitudinal crack were 
constructed to determine whether cracks generated by hypervelocity impact could 
cause catastrophic rupture due to internal pressure. The models gave results 
consistent with NASA-FLAGRO. Furthermore, equations were developed to 
allow linear elastic fracture mechanics analysis of a cylinder with cracks 
emanating from a hole, as well as influence (or weight) functions for arbitrary 
loads. Other finite element models were developed on CADDS to actually model a 
cylinder with stiffeners and a hole, but the time required to develop the model 
geometry and perform the minimum NASTRAN analysis prohibited completion 
of those analyses. 

Two different failure mechanisms, spall and necking, that were not originally 
anticipated, were predicted for an early Space Station Freedom wall design. 

These are discussed in the Handbook. Since the bumper and MLI plus beta cloth 
blanket were not optimized for the rear wall thickness, a higher loading rate on 
the rear wall is not well approximated by the Gausian momentum intensity 
distribution assumed by Wilkinson. This can lead to necking away from the 
centerline, or spall in the rear wall. Wide area spall of the rear wall can occur at 
velocities greater than 6 km/s and has been observed in 0.32 mm 2219 sheet. 
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3.3 Conclusions 

The measurements of rear wall thinning and deformation provide a valuable 
benchmark for transient elastic-plastic analysis codes. 

Thinning of the rear wall did not depend on edge constraints (although overall 
displacement was influenced by the edge constraints). 

The NASTRAN pie shaped sheet element did not produce results consistent 
with HULL or EPIC which were 2D axisymetric with 4 to 6 brick elements true the 
thickness. Some of the problem with the problem with the NASTRAN runs may 
be that bending was significant for the assumed momentum intensity distribution 
and rear wall thickness. (This was at the bounds of where Wilkinson noted the 
beginning of bending effects.) 

Three high velocity failure mechanisms are of interest: 

• General deformation, bulging and failure at the centerline as predicted by 
the Wilkinson’s model. 

• Plug or trapdoor type failure with necking and failure at the radius of the 
debris cloud. Heavier intermediate shields or catchers can slow down the 
momentum transfer to the rear wall transition to a Wilkinson failure 
mode. 

• Spall of a large area due to a impulse from the liquid or gaseous debris 
cloud. Spall failure should extrapolate with projectile energy rather than 
momentum. (The peak pressure varies linearly with impact velocity for 
the same momentum.) This is not favorable for PNP predictions, but 
should be beneficial to PNCD if more of the momentum passes through 
the wall. 

The linear elastic fracture mechanics stress intensity influence function for a 
longitudinal crack in a cylinder was calculated. It is very different than the 
influence function for a crack in a flat plate. Loads at the center of the crack have 
a significant bending effect not present in the flat plate. 

Crack growth due to impact may be analyzed by an energy balance between the 
initial kinetic energy of the rear wall, and energy dissipation by deformation and 
cracking. Estimating the kinetic energy of the rear wall is the most questionable 
part of the analysis. 

The critical flaw length that will continue to propagate in SSF (rear wall 
thickness = 0.125 inches (0.318 cm)) is calculated to be 23 cm in total length. While 
the central hole may make the critical crack length longer, due to air pressure on 
the inward-bent petals, this effect should not be counted on to prevent fracture. 

The SSF integrally machined stiffeners will not stop a crack which spans the 
diagonal between two nodes (53.2 cm center to center longitudinally along the 
cylinder). The stiffener spacing would have to be one third the current spacing 
(53.2 cm center to center) to be of value as a crack stopper. 
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TASK 4 - HYPERVELOCITY TESTS 

41 Objectives , | 

The objective of this task was to perform hypervelocity tests to verify the debris 
clouds and test an advanced shield concept and analyze interior penetration 
effects. Several specific tests were designed by MSFC and the raw test results 
were sent to the technical monitor at MSFC for analysis. 

4.2 Results 

The debris cloud results were documented and analyzed in Attachment 4. 
Attachment 8 contains test results that were sent to MSFC for analysis. 

4.3 Conclusions 

These tests added to a unique and valuable set of reference data on the post 
impact debris cloud. 

A textured bumper was not any more effective than a flat sheet in shattering a 
spherical impacting at 45° obliquity at 6.4 km/s. 

The interior penetration effects provided a unique capability for flash 
radiographic measurements of fragments in the interior of a module after an 
impact. 


TASK 5 - CONSULTATION TO MSFC 

5.1 Objectives 

The objective of this task was to provide consultation on related issues and 
problems as requested by NASA-MSFC. 

5.2 Results 

This included telecons, informal analyses and reports (included in 
Attachment 9), and the BUMPERII suggestion and problem reports (included in 
Attachment 3). 


TASK 6- HANDBOOK 

6.1 Objectives 

The objective of this task was to prepare a handbook for meteoroid and space 
debris damage assessment. 

6.2 Results 

The Handbook for Structural Damage Prediction and Analysis for 
Hypervelocity Impact is included as Attachment 7, 
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7 RECOMMENDATIONS 


7.1 Hypervelocity Testing 

Hypervelocity testing is expensive and restricted in velocity compared to design 
requirements, but the testing of one and two bumper shields continues to reveal 
the phenomenology of penetration, which allows the analyst to develop better 
predictive tools. Launcher development is also progressing with the promise of 
well characterized (if not well shaped) projectiles in the 10+ km/s regime. The 
following recommendations highlight those areas in special need of development. 

7.1.1 Scaled Velocity Testing 

Scaled velocity testing has been well developed to characterize the post impact 
debris cloud. The use of cadmium on cadmium or lead on lead have provided 
good benchmarks for hydrocode analysts. Further use of this approach is needed 
to model rear wall failure in the 10+ km/s regime. In particular, oblique impacts 
and the influence of intermediate layers (such as the Space Station Freedom 
thermal blankets) need to be characterized. 

While the cloud is in the correct physical state, the physics of rear wall impact 
also need to be correctly modeled. The possible rear wall failure mechanisms for 
liquid/vapor impingement must be clearly defined and modeled. The rate of 
application will influence whether spall can occur and the type of momentum 
failure and tearing. Schmitt et al observed that spall occurred for aluminum-on- 
aluminum impacts but not for cadmium-on-cadmium for the same 
dimensionless impact speed. However, this is still the most promising technique 
to qualify analytical models of bulging and tearing that may occur for oblique 
impacts at 10+ km/s. 

7.1.2 Multiple Bumper Penetration Resistance Testing 

Experimental design techniques are needed to improve the multi-shock 

technique. Whipple bumper systems have been developed through testing to 
understand the physical phenomena. However, the number of variables in a 
multi-shock bumper system precludes this type of approach. The variables 
include the materials, thickness and spacing of each layer. Two bumper systems 
have started to address these issues, but four-bumper systems have been 
restricted to identical evenly spaced bumpers. Experimental design techniques 
may be able to further optimize the design of multiple bumpers. 

Multiple bumpers are better at breaking up particles than single bumpers. 

But multiple bumpers require greater standoffs spread out momentum. 

7.1.3 Crack Growth Testing 

An accurate accounting of momentum is needed to develop a general 
predictive methodology for crack growth. This should include the use of 
momentum pendulums for tests at and beyond the threshold of penetration. 

Boundary conditions should be well characterized before drawing conclusions 
about crack growth from small specimens. 


8 



7.1.4 Round Robin Light Gas Gun Testing 

Bumper hole size data has indicated that technique differences can influence 
the results of light gas gun tests. A small round robin at different laboratories 
should be developed to quantify if there is an influence on rear wall penetration. 
While differences in the ballistic limit may be within normal scatterband, 
flattening of the projectile during launch may be a significant contributor to the 
scatterband. 

7.2 Impact Analysis 

7A2 Hydrocode Analysis 

Hydrocode predictions of rear wall momentum should be recorded both when a 
penetration is predicted and when it is not predicted. This is necessary for 
bulging and cracking analyses at much longer times than used for cratering or 
hole-out analyses. (Typical durations are 5 to 20 microseconds for penetration by 
cratering and 50-100 microseconds for bulging, and longer for tearing.) 

Large standoff multiple bumper analyses, both monolithic and cloth bumpers, 
present extreme challenges to hydrocode modelers. However, these analyses are 
needed for the next generation of bumper systems. 

7.2.2 Spall Analysis 

Rear wall spall typically needs finer meshes than used currently, but it is 
essential to account for this failure mechanism. Parametric studies of debris 
cloud interaction with the rear wall are needed to develop better models. Again, 
oblique impact is most significant. 

7.2.3 Deformation Analysis 

Continued deformation analysis is needed with newly developed post-impact 
cloud models and better material properties. The most success has been achieved 
with Gausian momentum intensity distributions as starting conditions, but a 
time dependent momentum transfer is more representative of the physics when 
an intermediate shield/catcher is present. 

Resolve the problem in using sheet elements for rear wall deformation (as 
available in NASTRAN and other codes). Examine the sensitivity of HULL and 
EPIC to the number of brick elements through the thickness. The simpler 3D pie 
slice model with sheet elements would be desirable for crack petaling studies 
instead of a 3D model with multiple brick elements through the thickness. 

7.2.4 Fracture Analysis 

R-curve analysis of crack growth should be developed to account for plastic 
behavior at the crack tip. The state of stress and strain around an impact induced 
crack is not the same as a fatigue crack. Elastic-plastic crack growth through an 
integrally machined stiffener is also different than the elastic solutions developed 
in this handbook. 
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7.3 Design Recommendations 

The use of crack stoppers should be used whenever practical. This includes 
add-on stiffeners as well as integrally machined stiffeners. Integrally machined 
stiffeners may be effective in containing the original impact induced flaw. Add-on 
crack stoppers could be selectively added to those areas at highest risk for 
momentum failure. Titanium crack stoppers between waffle nodes could pick up 
the load from a large petaled hole. 

7.4 Environment/Probability Analysis 

There are several improvements that could be made to either or both 
BUMPERII and SD.SURF. 

7.4.1 Collision Avoidance 

For critical damage estimates, the effect of collision avoidance on the effective 
debris flux should be included. 

7.4.2 BUMPERII 

Make separate rootfiles for Geometry and Response. Use the SuperTab 
filename as the root for the GEM file. 

Read input from a file (besides using Batch.com) (Included in Martin 
revisions to BUMPERII version 1.2a.) 

Maintain Macintosh compatibility in future updates. Include the 
suggested fixes (to avoid jumping into IF-THEN or DO loops) in the next 
release of BUMPERII. 

For Macintosh applications with either larger problems or reduced 
memory requirements, there should be an option to compile the 
RESPONSE, GEOMETRY and SHIELD modifies separately. 

7.4.2.1 BUMPERII - RESPONSE 

Add the ability to repeat a PID from memory 

Edit an existing PID. Select it by number and replace with another or 
recalculate the ballistic limit surface. 

Open an existing RESPONSE file as new, or append it to the current file 
for editing. 

Save a RESPONSE file, edit the data in memory and save as a different 
RESPONSE file. 

Add the CONTOUR feature of stepping through shield sizes. This would 
be useful with the “All” option in PJSURF. 

7.4.2.2 BUMPERII - GEOMETRY 

Add a graphic to the Macintosh version to show the rotation axes. 
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- Include a NASTRAN translator with BUMPERII. (Although the 
PATRAN translator in version 1.3 reduces the need for this option.) 

7.4.2.3 BUMPERII - SHIELD 

Add an option to average the diameter to penetrate (and the PNP) over a 
range of obliquities to account for granularity in the GEOMETRY model. 
(This will account for the fact that modules have curved surfaces rather 
than flat facets modeled in SUPERTAB or PATRAN.) This option should 
be selectable by either range or FID number so that truly flat surfaces will 
be modeled correctly. This would require input by the user of the typical 
angle subtended by the flat facets. 

Automatically report the breakdown of PNP, exposure area and shield 
mass among different PIDs that may be selected in one range. 

7.4.3 SDJSURF 

Smooth A_SURF and or P_SURF output to remove the effects of using 
large facets (15°) and large steps in threat velocity (45 threats in 180 
degrees). This will be a more realistic approximation of a curved surface 
such as a module. It will also aid in the interpretation of the results. 

7.4.3. 1 SD.SURF (FORTRAN version) 

Incorporate directly into BUMPERII. 

Analyze meteoroids separately, or both meteoroids and debris. 

Make one large SDJSURF program and automatically pass file 
information. 

Add contour and 3d charting to the Macintosh version. 


7.4.3.2 SDJSURF - AJ3URF 

Create new files with different ranges in the same session (read GEM file 
only once), or create multiple tables in one file (must also modify 
P_SURF). 

Add the text based contour map used in PJ3URF to AJ3URF. 

Generate the area arrays for a sphere automatically without 
GEOMETRY. (As done in AREA_Maker Macro.) The sphere is useful for 
a spacecraft that is not Earth oriented. 

Generate the area arrays as done in the AREAJMAKER EXCEL Macro. 
(This would not do the shadowing performed by BUMPERII. Adaptation 
of existing graphics packages to GEOMETRY may be more effective.) 
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7 .4.3.3 SD.SURF - P_SURF 

Put subtotals for each obliquity and each velocity (1 km/s increments) on 
the graph (and add same feature to PNP_Template). 

Include shield mass calculations as performed by SHIELD. 

Write PID descriptions to screen before asking which PID to use. 

Output only those PID descriptions used in an analysis. 

Process multiple selectable PIDs instead of just one or all. 
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